Statistics of spinons in the spin-liquid phase of CS2CUCI4 
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Motivated by a recent experiment on Cs2CuCi4^ we study the spin dynamics of the spin-liquid 
phase of the spin- 1/2 frustrated Heisenberg antiferromagnet on the anisotropic triangular lattice. 
There have been two different proposals^ for the spin- liquid phase of CS2CUCI4. These spin- liquid 
states support different statistics of spinons; the bosonic Sp(iV) large- N mean field theory 2 '* predicts 
bosonic spinons, while the SU(2) slave-boson mean field theory—*- leads to fermionic spinons. We 
compute the dynamical spin structure factor for both types of spin-liquid state at zero and finite 
temperatures. While at zero temperature both theories agree with experiment on a qualitative level, 
they show substantial differences in the temperature dependence of the dynamical spin structure 
factor. 
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The existence of a spin-liquid state in frustrated quan- 
tum magnets has been one of the central issues in the field 
of strongly correlated systems j&LSaLifliii Recent experi- 
mental studies of various frustrated magnetic compounds 
provided excellent opportunities to investigate the com- 
petition between geometric frustration and quantum fluc- 
tuations, and its role in the occurrence of quantum spin- 
liquid states— ^ Much interest in two-dimensional quan- 
tum spin-liquid states has evolved since Anderson pro- 
posed that the doped spin-liquid state may hold the 
key to the puzzles of high temperature superconductiv- 
ity in Cupratespl^ More recently the search for spin- 
liquid states has been extended to geometrically frus- 
trated quantum magnetsi 7 i 8 i 9 i 10 

The hallmark of the spin-liquid state is the exis- 
tence of fractionalized excitations. In this regard, 
the two-dimensional frustrated Heisenberg antiferromag- 
net CS2CUCI4 provides a useful realization of a two- 
dimensional spin-liquid state: A recent neutron scatter- 
ing experiment 1 showed the remarkable result that the 
dynamical spin structure factor <S(q, uS) does not exhibit 
well-defined peaks corresponding to spin-1 magnons. In- 
stead, there exists a continuum of excitations which has 
been interpreted as the indication of pairs of deconfincd 
spin-1/2 spinons in the underlying spin-liquid state. 

The minimal Hamiltonian describing CS2CUCI4 is ar- 
gued to be the spin-1/2 Heisenberg antiferromagnet on 
an anisotropic triangular latticed 
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where Si is the S = 1/2 spin operator on site i, and 
Ji, J2 are the exchange couplings along two different 
types of bonds, as indicated in Fig. 1. There have 
been at least two different proposals for the spin-liquid 
state of this model in connection to the experiment on 
CS2CUCI4. These proposals are based on two differ- 
ent mean field approaches: the bosonic Sp(7V) \arge-N 
mean field theory^ and the SU(2) slave-boson mean field 
theory— One of the distinguishing properties of the 
two resulting spin-liquid states is the statistics of the frac- 
tionalized excitations, i.e., spinons; the bosonic Sp(Af) 
large-TV mean field theory supports bosonic spinons 



while the SU(2) slave-boson mean field theory leads to 
fermionic spinons. 




FIG. 1: The antiferromagnetic exchange interactions in 
Cs 2 CuCl 4 with Ji = 0.125meV and J 2 = 0.375meVii 

The statistics of the spinons may be a useful character- 
istic of the underlying spin-liquid states. In particular, an 
important question is whether spin-liquid ground states 
with the same symmetry, but different spinon statistics, 
are necessarily distinct; this question has been only par- 
tially answered—^ In this paper, we would like to achieve 
a more moderate goal; while the above question about 
the ground state cannot be answered by investigating 
mean field theories, one may still be able to identify 
which mean-field state provides a more faithful approxi- 
mation of the true spin- liquid state of the material. If the 
mean-field state is a good representation, residual quan- 
tum fluctuations will be small, and are not expected to 
influence the system's responses in a qualitative way. We 
aim to achieve this goal by comparing the spin excita- 
tion spectra of the different proposals with the experi- 
mentally measured spin structure factor. We compute 
the dynamical spin structure factor for the spin-liquid 
states in both the bosonic Sp(iV) large- N and the SU(2) 
slave-boson mean field theories. We find that, at tem- 
peratures well below the temperature scale set by the 
exchange couplings, the temperature evolution of the dy- 
namical spin structure factor depends significantly on the 
spinon statistics: While the fermionic spinons are fairly 
insensitive to temperature, the bosonic spinon spectrum 
shows significant changes in this temperature range, as 
described below. 

A useful bosonic representation of the Hamiltonian in 
Eq. n can be obtained by generalizing the physical spin 
SU(2) = Sp(l) symmetry to Sp(iV). The generalized spin 
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operators can be expressed in terms of boson operators 
b\ a , ha on every site i, where a = 1 . . . 2N is a Sp(N) 
indexi^ The constraint b^ a b a = rib is imposed on each 
site to fix the number of bosons (rib = 25 for N = 1). In 
the mean- field theory, the Sp(iV) Hamiltonian is solved 
in the N — > oo limit with k = rib/N fixed A The mean- 
field phase diagram as a function of J2KJ1 + J2) and 
1/k contains both magnetic long-range-order (LRO) and 
short-range-order (SRO) phases (see Fig. 4 of Ref. Q) at 
zero temperature. For bare parameter values relevant to 
the material (k = 1 and J2/J1 — 3), the large- N phase 
diagram predicts a spin-ordered ground state with an in- 
commensurate wavevector in the J 2 -dircction. 2 However, 
finite-A r corrections will move the phase boundaries, so 
that the physical spin- 1/2 limit could in fact be described 
by the incommensurate SRO (spin-liquid) phase with de- 
confined spin- 1/2 spinons ji The material CS2CUCI4 ac- 
tually exhibits long-range order at temperatures below 
0.62K, which can be suppressed by an applied magnetic 
field. Since this ordering is due to the small interlayer 
coupling (J z < 10~ 2 J 2 )fi we expect the excitation spec- 
trum to be indistinguishable from the disordered state 
at energies above this (small) scale. In our strictly two- 
dimensional model, no ordered state can occur at any 
finite temperature. Here, we consider the point n = 0.64 
and J2/J1 = 3 in the large- N phase diagram (see Fig. 
4 of Ref. y), which lies in the SRO phase close to the 
LRO phase boundary, as the possible spin liquid state 
relevant to CS2CUCI4, which is experimentally known to 
be close to a quantum phase transition. The spinon dis- 
persion in this case has a small gap at q = (0.267T, 0.267r) 
of approximately 0.03J2, which is much smaller than the 
experimental resolution of about 0.5J 2 i 

On the other hand, spin-liquid phases can also be ob- 
tained from a fermionic representation via the SU(2) 
slave-boson approachi&£*i4 This approach utilizes a hid- 
den SU(2) gauge symmetry in the fermionic represen- 
tation of the Heisenberg models Introducing two SU(2) 
doublets, ip?! = (/it, /ij.) and ipf 2 = (f n , -f^), to rewrite 
the destruction operators of spin-up and down states^ 
the mean-field Hamiltonian can be expressed in terms of 
the 2x2 Hermitian matrices 

flmf = - ^2(i>huijipj a + h.c.) + aoV4 r Vfa- 

ij i 

Here, t 1 (I = 1, 2, 3) are the Pauli matrices, and the La- 
grange multiplier a l is used to enforce the constraint 

(V4a r V-ta!) = 0££ It has been proposed that the spin- 
liquid phase with the following mean-field ansatz is 
the most likely candidate for the spin-liquid phase of 
Cs 2 CuCl 4 ^ 

Ui, i+ x = X t2 i u i,i+y = -XT 2 , "M+z+y = At^ , 

a o' 2 = °> a o = fl 3- ( 2 ) 

Here the parameters x, A and 03 have to be determined 
self-consistently. According to the classification scheme 



of Ref. 0, this is one of the possible U(l) spin-liquid 
phases on the anisotropic triangular lattice. The exci- 
tation spectrum has gapless points^ at (0, 0), (71", n), and 
(q, q), where q — ±ir/2±e, with e = 0.047T. Although the 
mean-field solution is expected to be modified by gauge 
field fluctuations, it may still provide a qualitatively cor- 
rect description of the spin excitation spectrum. 
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FIG. 2: Intensity plots of x"(q, tJ )> with q oriented along the 
direction of J2, computed in (a) the Sp(iV) and (b) the SU(2) 
approach. Darker areas indicate higher scattering intensity. 
Top panels represent the results at T = and bottom panels 
indicate the results at T — O.I5J2. 



The spin-liquid states obtained from the two ap- 
proaches differ in important aspects: one of them is a Zi 
spin-liquid phase with gapped bosonic spinons (although 
the gap is small), and the other is a U(l) spin-liquid 
state with gapless fermionic spinons. It is therefore in- 
teresting to ask which spin liquid state is more likely to 
describe CS2CUCI4. We therefore calculate in both cases 
the dynamical spin structure factor defined as 5(q, w) = 
-i(l + n(w))x"(q,w), where n(ui) = \ j(e^l T - 1) is the 
Bose thermal factor, and x"(q? w ) is the imaginary part 
of the dynamical spin susceptibility: 



X "(ci,oj)=3Jm\-J dte^e(t)([S z (<i,t),S z (^0)})) , 

(3) 

where spin rotation invariance has been used to write 
X = %Xzz- 
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In the Sp(N) approach x"(q, oj) is obtained as 



v" 

ASp(JV) 



= 3E 



(i + A(q))(n(Bk) - n(££))<J(w + £ k - £ k ) 



+(- + i4(q))(n(^) - n(£ k ))c5(u, + £ k - £ k ) 
+ (i - A(q))(l + n(£ k ) + n(KMu> + £ k + £ k ) 



(--4(q))(l + n(4)+n(£[))i( W 



(4) 



where A(q) = = C' 



and A! = 



A k + q . Here n(Ek) is the Bose thermal factor, £? k = 
A 2 — A 2 , is the bosonic spinon dispersion^, and A k = 
JiQi(smk x + sinfcj,) + J2Q2 sin^ + k y ). Qi and Q2 
are mean-field bond variables for nearest-neighbor and 
next-nearest-neighbor bonds, respectively, and A is a La- 
grange multiplier enforcing the constraint on the number 
of bosons. The temperature dependent values of Q%, Qi-, 
and A are determined self-consistently. 

In the SU(2) slave-boson approach, x" is given by 



XSU(2) 



(i + B(q))(/(ftO - /(££))*(" + £ k - £i) 



+(- + B(q))(/(££) - f(fr))8(w + £i- ftj 
-(i - S(q))(l - /(£ k ) - f(£L))S(cu + £ k + 4) 
+ (i - S(q))(l - /(ft k ) - f{£i))5(oj - ft, - E$ 



(5) 



k+q, and £> k 



where B(q) = k , £ k = ft 

Here, ft, = 2^/e 2 + D 2 is the dispersion of the fermionic 
spinons, 3 e k = A cos{k x +k y )— 03, D k = x(cos k x — cos 
and f{£k) is the Fermi distribution function. We assume 
that the values of the order parameter fields x an d A, 
which were obtained in Ref.y, do not vary significantly in 
the temperature range considered here, and we explicitly 
verified that our results are insensitive to any changes in 
03 with temperature. 

Our results for x"(q, a>) for the two spin-liquid states 
are shown in Figs. 2-4. The spin structure factor <S(q, u>) 
differs from x"(q, a;) by an overall thermal factor (1 + 
n(u>)). As a result the low energy spectral weight in 
x"(q, w) is somewhat suppressed compared to «S(q, oj), 
but the overall intensity distribution looks very simi- 
lar. To make direct contact with experiment, we show 
x"(q, oj) for both spin liquid states in Fig. 2, with q ori- 
ented along the direction of Ji (see Fig. 1). Figures 3 
and 4 show the same quantity along four scan trajecto- 
ries that were used in the experiment (see Figs. 2 and 3 
of Ref. [3). Note that the point k = 2n/b in Fig. 2a of 
Ref. 0] corresponds to q — n in Fig. 2 here. 
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FIG. 3: Scans of x"(q, ^) at various temperatures in the 
bosonic Sp(iV) approach along the scan directions A, B, C 
and D as used in Fig. 2(a) of Ref. The results at T = 0, 
O.I5J2 and O.2J2 are shown by the solid, dotted, and dashed 
lines, respectively. 
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FIG. 4: Scans of x"(q, w) at various temperatures in the 
fermionic SU(2) approach along the same four scan directions 
as in Fig. 3. The results at T — 0, O.I5J2 and O.2J2 are shown 
by the solid, dotted, and dashed lines, respectively. 



Both spin-liquid states show a continuum of spin exci- 
tations, which is a hallmark of spin- liquids with spin-1/2 
spinons. Moreover, both theories indicate strong scatter- 
ing around (7r/2, 7r/2) in agreement with the experiment. 
At zero temperature, the lower edge of the excitation 
spectrum in both theories has minima at (0,0), (tt, tt) 
and close to (7r/2,7r/2). The spinon spectrum in the 
Sp(N) approach has a small gap (0.03 J2) at the single in- 
commensurate minimum while it has two gapless incom- 
mensurate points in the SU(2) theory. These slight dif- 
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ferences, however, are below the experimental resolution 
(about 0.5J2)i The upper boundary of the scattering 
spectrum in the SU(2) approach is closer to the experi- 
mental results. On the other hand, the Sp(A) approach 
describes some aspects of the experimentally measured 
spectra at lower energy (see Fig. 2a of Rcf. [l|) better: 
there is much less scattering intensity around (tt, tt) in 
the Sp(iV) approach compared to the results of the SU(2) 
theory, and the peak intensity is close to the lower edge 
as opposed to the upper edge in the case of the SU(2) 
approach. The details of the spectra described above are 
expected to change, however, once fluctuations about the 
mean-held solutions are included. In addition, it is dif- 
ficult to distinguish the two spin-liquid states from their 
low energy excitation spectra under current experimen- 
tal resolution. Therefore, it is fair to say that, at zero 
temperature, both theories agree reasonably well with 
experiment on a qualitative level. 

At finite temperatures, however, one expects that the 
different spinon statistics will give rise to substantial 
differences in the spin excitation spectrum. Indeed, as 
shown in Figs. 2-4, our results at finite temperatures in- 
dicate that the scattering intensity extends to a broader 
range in the phase space of q, w, and \'' '(Qi w ) shows very 
different behavior with increasing temperature in the two 
approaches. This difference can be most clearly seen by 
comparing the temperature evolution of the intensities 
along the four scan directions used in the experiment 
(see Fig. 2 of Ref.Q). As shown in Fig. 3, the maximum 
intensity increases in the Sp(N) approach, and the exci- 
tation spectrum becomes more narrowly peaked around 
the maxima as temperature increases, while the overall 



intensity remains roughly constant. Additionally, a two- 
peak structure develops, which resembles the lineshape 
predicted by a spin-wave calculation in the LRO phase 
(see Fig. 3 of Ref.yj), and is presumably due to the vicin- 
ity of an ordered quantum ground state. On the other 
hand, as shown in Fig. 4, the evolution of the spin struc- 
ture factor with temperature is much less significant in 
the SU(2) theory, where only a very slight shift of spectral 
weight to lower energies is observed. Since the different 
temperature dependence arises from distinct statistics of 
the spinons, these results will not be affected fundamen- 
tally by fluctuations about the mean-field ground states. 

In summary, we studied the dynamic spin structure 
factor of different spin-liquid states obtained from the 
bosonic large- iV Sp(iV) and the fermionic SU(2) mean- 
field theories of a two-dimensional frustrated spin-1/2 
Heisenberg antiferromagnet, as a model for Cs 2 CuCi4. 
We found that at zero temperature the dynamic spin 
structure factor of both spin-liquid states compares fa- 
vorably with experiment. At finite temperatures, fun- 
damentally different behavior arises due to the different 
spinon statistics. The signatures we found in the tem- 
perature dependence of the spin structure factor can be 
used to compare both theories with future experiments, 
in order to determine which theory is best suited to de- 
scribe the spin liquid state of CS2CUCI4. More details 
will be presented in a future publication. 
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